Predicting decoherence in discrete models 
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The general aim of this paper is to supply a method to decide whether a discrete system decoheres 
or not, and under what conditions decoherence occurs, with no need of appealing to computer 
simulations to obtain the time evolution of the reduced state. In particular, a lemma is presented 
as the core of the method. 



I. INTRODUCTION 

At present, the study of quantum decoherence has acquired a central position in the theoretical research on quantum 
£NJ • mechanics. Although the orthodox environment induced decoherence (EID) approach addresses decoherence in open 
systems ((H-0)) many authors have stressed that close d sy stems may also experience the phenomenon of decoherence 
^ ' (13" Q3; we have worked from this perspective in [l8|-|29l|). 

In order to show that the two approaches must not be conceived as rival or alternative, but rather as complementary, 
we have developed a general theoretical framework for decoherence ([3(J, [3l|, [IH), which encompasses decoherence 
in open and closed systems. According to this general framework, decoherence is just a particular case of the general 
phenomenon of irreversibility in quantum mechanics (jH|, [HI)- I n fact, since the quantum state p(t) follows a unitary 
evolution, it cannot reach a final equilibrium state for t — > oo. Therefore, if we want to explain the emergence of 
non-unitary irreversible evolutions, a further element has to be added: we must split the whole space O of all possible 
observables into a relevant subspace Or C O and an irrelevant subspace. Once the essential role played by the 
selection of the relevant observables is clearly understood, the phenomenon of decoherence can be explained in three 
general steps: 

3 . 

Q" 1 1. First step: The space Or of relevant observables is defined 

2. Second step: The expectation value (Or)^, for any Or £ Or, is obtained. This step can be formulated in 
two different but equivalent ways: 

^ ! 

• (OR,) p (t) is computed as the expectation value of Or in the unitarily evolving state p(t). 
fNl ■ • A coarse-grained state PR(t) is defined by 

m 

O- (0 R ) p(t ) = {0 R ) PR(t) VO R eO R (l) 

' and its non-unitary evolution (governed by a master equation) is solved. 

3. Third step: It is proved that (Or) p ( t ) = (0 R ) PR ( t ) reaches a final equilibrium value: 



in 



- 1—^ 

X 



lim (0 R ) p(t) = lim (0 R ) PR{t) = (Or) p , = (Or) Pr , VOr € Or (2) 

The final equilibrium state is obviously diagonal in its own eigenbasis, which turns out to be the final equilibrium 
decoherence basis. But, from eq. ([2]) it cannot be concluded that lim^oo p(i) — p„; the mathematicians say that the 
unitarily evolving quantum state p(t) of the whole system has only a weak limit, symbolized as: 

W - lim p(t) = p» (3) 

t— >oo 

and equivalent to eq. @. Physically this weak limit means that, although the off-diagonal terms of p(t) never vanish 
through the unitary evolution, the system decoheres from an observational point of view, that is, from the viewpoint 
given by any relevant observable Or £ Or. From this general perspective, the phenomenon of decoherence is relative 
because the off-diagonal terms of p(t) vanish only from the viewpoint of the relevant observables Or € Or. 

This general framework strictly applies when the limit of eq. (J5|) exists. This happen when the off-diagonal terms 
of the density matrix vanish by destructive interference according to the Riemann-Lebesgue theorem. In fact, when 
the relevant observables Or G Or read 

/>oo />oo f. oo 

Or = \ 0(w)\uj) + / / 0(u>,L>')\uj,Lj')du)dL>' (4) 
Jo Jo Jo 
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where {|w}} is the eigenbasis of the Hamiltonian, \ui) = \uj)(uj\, \uj,uj') — and {|w), \uj,lu')} is a basis of Or, 

and the states are linear functionals belonging to 0' R , the dual of Or, 

/•OO />oo />oo 

p = duj + / p{uj,uj )(u),u) | djjjdjjj' (5) 

Jo Jo Jo 

where {(w|, (w,c</| } is the basis of 0' R , then the expectation value of any observable Or G O/j in the state pit) can 
be computed as 

poo /* oo />oo , 

(Or) P[t ) = W)0{lo) du + / p(u,u»)0(u,u') e 1 ^* dudu' (6) 
Jo Jo Jo 

When the function p(w, lu')0{uj, u/) is regular (precisely, when it is Li in variable v = lo —oj'), the Riemann-Lebesgue 
theorem can be applied to eq. ©. As a consequence, the second term vanishes and (0_r} p (() converges to a stable 
value 

poo 

(Or) pW — > / ^)0(.) = (O fl ) p , (7) 
Jo 

where is diagonal in the eigenbasis of the Hamiltonian. 

It is clear that the Riemann-Lebesgue theorem strictly applies only in cases of continuous energy spectrum. However, 
we also know that we can use the results coming from the continuous realm in quasi-continuous cases, that is, in 
discrete models where (i) the energy spectrum is quasi-continuous, i.e., has a small discrete energy spacing, and (ii) 
the functions of energy used in the formalism are such that the sums in which they are involved can be approximated 
by Riemann integrals. These conditions are rather weak: in fact, the overwhelming majority of the physical models 
studied in the literature on dynamics, thermodynamics, quantum mechanics and quantum field theory are quasi- 
continuous. 

The general aim of the present work is to supply a rigorous formulation of this intuitive idea. In particular, we will 
develop a discrete analogue of the Riemann-Lebesgue theorem, and this task will lead us to introduce a lemma in terms 
of which it is possible to predict whether a discrete system decoheres or not. For this purpose, the paper is organized 
as follows. In Section 2 we will consider the three cases that can be distinguished regarding the discrete analogue 
of the Riemann integral involved in the Riemann-Lebesgue theorem. Section 3 will be devoted to formulate the 
conditions for the validity of the discrete analogue of the Riemann integral. In Section 4, the lemma that constitutes 
the core of the proposed method is presented, and its relationship with the Discrete Fourier Transform is pointed 
out. Finally, in the Conclusions we will stress the fact that, whereas the usual strategy in this field is to rely on that 
computer simulations, our lemma makes possible to draw conclusions without that strategy; moreover, our results 
are particularly adequate to take advantage of the many mathematical methods of software engineering based on the 
Discrete Fourier Transform. 



II. THREE CASES IN THE DISCRETE ANALOGUE 



The Riemann-Lebesgue theorem —the mathematical expression of destructive interference— establishes that 

f{v) E h 1 lim / dvf{v)e ivt = (8) 

In order to obtain a version that can be used in the discrete case, let us analyze the discrete analogue of the Riemann 
integral R(t): 

R(t) = ^ duf(u)e^ — > R D (t) = J2 jjf (jj) ^ (9) 

where < j/N < 1, and t is a dimensionless time. Since the function Ro(t) is a finite sum of sine functions / (4r) e*»', 
it has a recurrence or Poincare time. Given an initial state Rd(0), the Poincare time tp is defined by i?_o(0) = Roitp). 
This means that 

N 1 / ' \ 

£ N f (jf) tp = ^ (10) 

Since the system comes back to the initial state when t — tp , there is no rigorous discrete analogue of the 
Riemann-Lebesgue theorem. Nevertheless, three possible situations can be distinguished: 
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1. If N oo, then — ^| becomes infinitesimal and tp — > oo. Therefore, this situation can be considered a 
continuous-spectrum case where the Riemann-Lebesgue theorem can be applied. 

2. If N is large, then | — | is very small. Then, the sum turns out to be close to the Riemann integral, and the 
the situation can be approximated to the continuous-spectrum case where the Riemann-Lebesgue theorem can 
be applied. This condition is satisfied in a concrete example in [35j : in spite of the fact that, strictly speaking, 
a system with discrete spectrum never reaches equilibrium due to Poincare recurrence, that paper shows that, 
for times t tp, the discrete spectrum can be approximated by a continuous spectrum where the considered 
functions satisfy the usual conditions of regularity and integrability. 

3. If N is not large, then | — | is far from being infinitesimal, and the sum cannot be approximated by a 
Riemann integral. Consequently, the Riemann-Lebesgue theorem is not applicable because there is no destructive 
interference. 

Let us notice that the difference between cases 2 and 3 is not absolute, to the extent that the precise criterion to 
decide when N is large has not been defined. In the following subsections such a criterion will be established. 



III. CONDITIONS FOR THE VALIDITY OF THE DISCRETE ANALOGUE 



The problem is to find the conditions for the time tp such that Rd^) when t — > tp: therefore, in the time-scale 
[0,£f] it can be considered that N is large enough to make the continuous-spectrum approximation applicable. In 
order to face this problem, we begin by consider a fixed time t and by decomposing the exponential of eq. Q as 



N N „ N 



Rd ® = 12 if (*<) elXtt = E if (*<) cos (***) + 1 E if (*«) sin o***) ( n ) 



i=0 i=0 i=0 



where the points Xi = i/N belong to a discrete set {x{\ with i = to N. Let us analyze the particular case where 
f(Xi) = l, precisely, 



N -, N 



r{ d (*) - E i cos ^ + 1 E i sin ^ ( 12 ) 



N " ' ^— ' N 

i=0 i=0 



In particular, we will consider the sums 

N N 

-^cos(^) = i4 1C V) -^sin^t) = < 1S) (t) (13) 



N ^ v ' u w N 

i=0 i=0 



We begin by considering w D '(t), because the case of R D (t) will be analogous. 

The sum R D (t) vanishes when its terms cancel by pairs, that is, when for any Xi € {x{\ 

cos (xjt) + cos (xit + 7r) = (14) 

where Xk — X4 + ir/t G {xi}. Since Xi — i/N , this happens when 

k i 7T , 7T . 

— = — + - => k - i = -N (15) 

N N t t V ' 

where k — i e N and N 6 N. However, since in general ir/t ^ N, the condition of eq. (fTB"]) is not always satisfied. 
Then, instead of requiring that the terms of R D (t) cancel with each other exactly, we will only require that the 
corresponding difference be small in the following sense: 

|cos (xit) + cos (xit + 7r + Sj)\ < e <C 1 (16) 

where now x^ = Xi + n/t + Sj/t G {x^ and 

Sj = min {Si} with Si — x^t — x{t — n (i — 0...N) (17) 
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If 5j < 1, the Taylor development of cos (x{t + tt + Sj) leads to 

|cos (xit) + cos (xit + tt + Sj)\ ~ | sin (xtf) Sj + cos (x*t) 8A<e (18) 

But, on the other hand, 

| sin (x{t) Sj + cos (xit) Sj \ = \Sj \ |sin (xit) + cos (xit) Sj \ 

< \Sj\ (|sin(x i t)| + IcosOji)! \Sj\) < \Sj\ (1 + \5j\) < \Sj\ (19) 

Then, if \5j\ < e -C 1, from eqs. (IT8"|) and (TT9")) we obtain the condition of eq. (IT5|) . Therefore, the condition of 
approximate cancelling is (see eq. (|17p) 

i5j| < e <C 1 with Sj = Xkt — Xjt — tt (20) 

Now we will express the condition of approximate cancelling of eq. (1201) in terms of the time t. Let us begin by 
noticing that the condition is not satisfied for t — 0, since t = Sj = —tt =>■ \Sj\ > e. Then, the first condition is 
t > 0. Now, by recalling that X{ = i/N, from the expression of Sj in eq. ([20]) we obtain 

SjN 

k- 3 = — + ^ r (21) 



But since j, k £ N, and j, k £ [0, N], then for j < k 

l<k- 

Then, if \Sj \ < e <T 1, eq. implies that 



l<fc-j<AT =► 1<^ + ^<7V (22) 



7rA^ 

1 < —j— < N (23) 



Therefore, the condition \5j - \ < e <C 1 of approximate cancelling turns out to be 

tt < t < ttN (24) 

Up to this point we have proved that, for tt < t < ttN, cos (xjt) approximately cancels with cos (x j t + tt + Sj ) . 
Figure ?? shows an example of this situation, where point 1 is cancelled by point 9, point 2 by point 10, ... , point 8 
by point 16. However, this is not the most general case, since the points cancel by pairs only when t = Ittu. In the 
general case there are points with no counterpart to be cancelled. An example of this situation is shown in Figure 
??, where points 13, 14, 15 and 16 are not cancelled. 

In order to analyze this general situation, let us consider the "worst" case, when the points of a whole half-period 
are not cancelled, namely, when t = (2n + 1) tt. Since in t there are N + 1 points, in such a half-period there are 
(N + 1) /(2n + 1) points, whose contribution r n (t) to the sum R D (t) is 

N+l 1 
^ 2n + l 

r *^=N ^ c °s(^) (25) 

The upper boundary of this contribution is 

N + l 1 

1 2 , s 1 N+l 1 TT 

TV ^2/1+1 2n + 1 t 

Then, the contribution rv(t) of the non cancelled points is irrelevant when r w (t) = ir/t < e <C 1, and this adds the 
condition 

t > tt (27) 

Summing up, if we combine eqs. (1241) and (|27[) , we obtain the condition on the time-scale that guarantees the 
approximate cancelling of the terms of Rp C ^ (t) : 

If tt < f < ttN => R { D C \t) < e < 1 (28) 

By means of the same argument applied to Rp S \t) we obtain an analogous result that, when combined with eq. (|28|). 
leads to 

If tt < f < ttN R ( ^(t) < e < 1 (29) 
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FIG. 1: Point 1 is cancelled by point 9, point 2 by point 10, ... , point 8 by point 16. 

IV. A LEMMA FOR THE APPLICATION OF THE DISCRETE ANALOGUE 

Up to this point we have studied the case / (xj) — 1; now we will consider the case / (x^ ^ 1. In order to compute 
Po(i) as defined by eq. (fTT|) . we have to ask a certain degree of regularity to function f {xi). The first step consists 
in splitting the set {xt} in G subsets of (P + 1) consecutive points: 

G G 

i x i) = IJ {x(k-i)(p+i)+i, -,%k{P+i)} = U X k (30) 
fe=i fc=i 

Now we relabel the points Xj G X k : since j — (k — 1) (P + 1) + 1 to k(P + 1), we can replace the index j by the 
index = j + (1 — fc) (P + 1) — 1, and we obtain 

Xj e X k — > x rk e X k with r k = 0, P (31) 

Then, we define 

Definition 1 . Let {xi} be a set of points uniformly distributed (or equidistant), with i € [0,iV] and N 1. 
The set {xi} is said to be quasi- continuous of class 1 if 3G £ N, BP G N such that P 1 and {x;} = 

Ufe=i { a -(fe-i)(P+i)+ij •••) x k(P+i) } — Ufc=i-^fc- ^ e se ^ ^ * s called the k component of the quasi- continuous de- 
composition. 

If the function f(x rk ) is almost constant in Xk, i.e. 

f(x rk ) = C k (32) 

Then, we can define 

Definition 2 . Let f(xi) : R — > M 6e a discrete function defined over the quasi- continuous set {x{\ of class 1. If for 
every component X k of a quasi- continuous decomposition f{x rk ) = C k , with x Tk G X k , we say that /(x,) G C\ . 
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FIG. 2: Points 13, 14, 15 and 16 are not cancelled. 



Therefore, when f{x{) G C\, the discrete function can be written as 



*° w - E (*o eiXit = E | f E ) = E ^ (E ^ 

i=0 fc=l \r fc =0 / fe=l \rs.=0 



If we define the function 



4 fc) (*) = E^ 



e k 



then the discrete function Ro{t) results 



r-fc=0 



P 



fc=i 



(33) 



(34) 



(35) 



Under this form, the condition of eq. (|29[) obtained in the previous subsection can be applied to each Rp (t) : 

If 7T < t < ttP => R ( D(t) < e < 1 (36) 
When this condition is satisfied, the sum i?£>(£) results 



p p p pr^ 

= E v^ffw < E ^ E v Ce = ^r Ce = Ce 



(37) 



fc=i fc=i fc=i 

where e = maxfe {sk} and C = maxt {Cfc}. As a consequence, if we consider that tp — 27r, we have proved that 
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Lemma 1. Let f{xi) be defined over the quasi- continuous set {x^} of class 1, with i — 1...N. If f(xi) s C\, then 

- 1 

There arc different kinds of functions for which the sum Ro{t) vanishes and that could be characterized by further 
lemmas, but we will not consider those cases now. Nevertheless, a practically useful remark is in order. The condition 
of eq. (|32l) (that the function f(xi) be approximately constant in each element Xk of the partition) can be expressed 
under a matematically more elegant form. Given a function f(xi), its Discrete Fourier Transform (DFT), as used in 
signal analysis ([36j]-[39j]), is defined as 

N 1 

/(<)=Eiv /(a;i)efa ' t (39) 

i=0 

This result may be very useful in practice, in particular in cases in which Lemma 1 is difficult to be applied. In fact, 
when we realize that the sum Ro(t) corresponding to the function / (xi) is precisely the DFT of / (xi), we can use 
all the properties of DFT —as linearity, symmetry, time-shifting, frequency-shifting, the time-convolution theorem, 
and the frequency-convolution theorem— to study Ro(t). Moreover, we can take advantage of the large amount of 
software designed to compute DFT and profusely used in physics and engineering. All these resources, which are 
standard tools in signal analysis, may prove to be extremely useful for studying decoherence in discrete models. 

V. CONCLUSIONS 

In this paper we have offered a discrete analogue of the Riemann-Lebesgue theorem and, on this basis, we have 
introduced a lemma relevant for discrete models, which provides a criterion for deciding whether or not the system 
decoheres with no need of numerical simulations. Moreover, we have shown how the large amount of mathematical 
methods of software engineering based on the Discrete Fourier Transform can be used to predict decoherence in 
discrete models. 



VI. ACKNOWLEDGMENTS 

This research was partially supported by grants of the University of Buenos Aires, the CONICET and the FONCYT 
of Argentina. 



[1] H. D. Zen, Found. Phys., 1, 69, 1970 

[2] H. D. Zeh, Found. Phys., 3, 109, 1973. 

[3] W. H. Zurek, Phys. Rev. D, 26, 1862, 1982. 

[4] W. H. Zurek, Progr. Theor. Phys., 89, 281, 1993. 

[5] J. P. Paz and W. Z., "Environment-induced decoherence and the transition from quantum to classical", in Dieter Heiss 
(ed.), Lecture Notes in Physics, Vol. 587, Springer, Heidelberg-Berlin, 2002. 

[6] W. H. Zurek, Rev. Mod. Phys., 75, 715, 2003. 

[7] L. Diosi, Phys. Lett. A, 120, 377, 1987. 

[8] L. Diosi, Phys. Rev. A, 40, 1165, 1989. 

[9] G. J. Milburn, Phys. Rev. A, 44, 5401, 1991. 
[10] R. Penrose, Shadows of the Mind, Oxford University Press, Oxford, 1995. 
[11] G. Casati and B. Chirikov, Phys. Rev. Lett, 75, 349, 1995. 
[12] G. Casati and B. Chirikov, Physica D, 86, 220, 1995. 

[13] S. Adler, Quantum Theory as an Emergent Phenomenon, Cambridge University Press, Cambridge, 2004. 
[14] R. Bonifacio, S. Olivares, P. Tombesi and D. Vitali, Phys. Rev. A, 61, 053802, 2000. 
[15] G. W. Ford and R. F. OGonnell, Phys. Lett. A, 286, 87, 2001. 
[16] M. Frasca, Phys. Lett. A, 308, 135, 2003. 

[17] A. C. Sicardi Shifino, G. Abal, R. Siri, A. Romanelli and R. Donangelo, "Intrinsic decoherence and irreversibility in a 
quasiperiodic kicked rotor", arXiv:quant-ph/0308162 2003. 



8 



[39] 



M. Castagnino, Int. J. Theor. Phys., 38, 1333, 1999. 

M. Castagnino and R. Laura, Phys. Rev. A, 62, 022107, 2000. 

M. Castagnino and R. Laura, Int. Jour. Theor. Phys., 39, 1767, 2000. 

M. Castagnino and O. Lombardi, Int. Jour. Theor. Phys., 42, 1281, 2003. 

M. Castagnino, Physica A, 335, 511, 2004. 

M. Castagnino and A. Ordonez, Int. Jour. Theor. Phys., 43, 695, 2004. 

M. Castagnino, Braz. Jour. Phys., 35, 375, 2005. 

M. Castagnino and O. Lombardi, Phys. Rev. A, 72, 012102, 2005. 

M. Castagnino and O. Lombardi, Phil. Scie., 72, 764, 2005. 

M. Castagnino, Phys. Lett. A, 357, 97, 2006. 

M. Castagnino and M. Gadella, Found Phys., 36, 920, 2006. 

M. Castagnino and O. Lombardi, Chaos, Sol. Frac., 28, 879, 2006. 

M. Castagnino, S. Fortin, R. Laura and O. Lombardi, Class. Quantum Grav., 25, 154002, 2008. 
M. Castagnino, S. Fortin and O. Lombardi, Jour. Phys. A: Math, and Theor., 43, 065304, 2010. 
M. Castagnino, S. Fortin and O. Lombardi, Mod. Physics Lett. A, 25, 1431, 2010. 
R. Omnes, "Decoherence: an irreversible process", arXiv:quant-ph/0106006 t 2001. 
R. Omnes, Phys. Rev. A, 65, 052119, 2002. 

F. Gaioli, E. Garcia- Alvarez and J. Guevara, Int. J. Theor. Phys., 36, 2167, 1997. 

S. Bochner and K. Chandrasekharan, Fourier Transforms, Princeton University Press, London, 1949. 

E. Oran Brigham, The Fast Fourier Transform and its Applications, Prentice-Hall, New Jersey, 1988. 

W. Press , J. Teukolsky , E. Vetterling and B. Flanery, Numerical Recipes in Fortran, Cambridge University Press, 

Cambridge, 1994. 

J. Kauppinen and J. Partanen, Fourier Transforms in Spectroscopy, Wiley- VCH, Berlin, 2001. 



